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Abstract: Many models that include small extra space dimensions predict graviton 
states which are well separated in mass, and which can be detected as resonances in 
collider experiments. It has been shown that the ATLAS detector at the Large Hadron 
Collider can identify such narrow states up to a mass of 2080 GeV in the decay mode 
G —> e + e _ , using a conservative model. This work extends the study of the e + e - channel 
over the full accessible parameter space, and shows that the reach could extend as high as 
3.5 TeV. It then discusses ways in which the expected universal coupling of the resonance 
can be confirmed using other decay modes. In particular, the mode G —> 77 is shown 
to be measurable with good precision, which would provide powerful confirmation of the 
graviton hypothesis. The decays G —> R + R~, W + W~, Z°Z° and jet-jet are measurable 
over a more limited range of couplings and masses. Using information from mass and 
cross-section measurements, the underlying parameters can be extracted. In one test 
model, the size of the extra dimension can be determined to a precision in length of 
7 x 10 ~ 33 m. 
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1. Introduction 

An exciting idea to be tested in high-energy collider experiments is the possible existence 
of narrow graviton resonances in the TeV energy range. Such resonances are predicted in 
models with small extra spatial dimensions. An example is the localized gravity model of 
Randall and Sundrum (RS) pj. This model aroused great theoretical interest because it 
motivates the weak-Planck scale hierarchy via an exponentially suppressed warp factor in a 
non-factorisable geometry. Many possible extensions and elaborations of this type of theory 
are being discussed in the literature 1-0Q-0B0B1-0- Problems with negative tension 
brane instability in the original RS model are solved in some other models. For example, 
placing the branes on fixed points of orbifolds projects out the negative energy modes [|TTfl. 


- 1 - 





























An extra scalar with couplings on the branes can be used to naturally stabilise the brane 
separation |2j]. Warped extra dimensions (with associated graviton resonances) have also 
been considered in the context of supersymmetry m 0 m [l6| . Thus the discovery of 
TeV-scale graviton resonances remains a possibility that needs to be considered seriously 
in preparing for future collider experiments. 

In [17], the detection of a narrow graviton resonance using the ATLAS detector at the 
Large Hadron Collider was considered. The main aim of that paper was to establish the 
discovery limit in the most favourable decay channel, G —> e + e”. A similar study for the 
CMS detector has since been reported 


In |17| the angular distribution of the lepton 
pair was also studied, and it was shown that the spin-2 nature of the resonance could be 
confirmed up to a mass somewhat below the discovery limit. 

Apart from its unique spin, the most striking characteristic of the graviton is its uni¬ 
versal coupling to all types of matter and gauge fields. In the present paper we consider the 
accuracy with which the couplings of a narrow graviton resonance to leptons, electroweak 


bosons, hadronic jets and Higgs bosons could be measured at the LHC. As in [17], we use 
the expected properties of the ATLAS detector as a guide to experimental limitations and 
the simplest RS model to characterise the resonance parameters, but our results should 
apply to other general-purpose detectors and to a broad class of models. We do, however, 
assume that all matter and gauge fields are confined to the physical brane and do not 
propagate into extra dimensions, thus excluding models of the type considered in Jlf], 20]. 


In the simple RS scenario, a 5-dimensional non-factorizable geometry is used, with two 
3-branes of opposite tension. A graviton Kaluza-Klein spectrum is created, with a scale 


At r — Mpie 


—kr c TT 


( 1 . 1 ) 


where Mpi is the reduced effective 4-D Planck scale, r c is the compactification radius of 
the extra dimension, and k is a scale of the order of the Planck scale. The geometrical 
exponential factor (the ‘warp factor’) generates TeV scales from fundamental Planck scales 
and hence motivates the weak-Planck hierarchy, if kr c ~ 12. 

The masses of the graviton resonances are given by 


m r 


= kx r , 


—kr c 7r 


= x n {k/Mpi) A„ 


( 1 . 2 ) 


where x n are the roots of the Bessel function of order 1 (x n = 3.8317, 7.0156,10.1735 for 
n = 1,2,3). The massive graviton excitations couple with equal strength to the visible 
sector ]^l]. However, the higher modes being suppressed by the falling parton distribution 
functions, only the lightest mode is considered in this paper. This does not in any way 
affect the generality of the approach, as the analysis can be applied to any such resonances, 
including the higher modes, so long as the resonances are narrow and sufficiently separated 
from the other modes. This is in contrast to studies in which many excitations, each with 


small coupling, contribute to some scattering process |2^, |2^, |24|, |25|, |(| . For brevity, we 
refer to the first massive resonance, with mass me = mi, as “the graviton”. 

In the RS model, the couplings of the graviton are given by 1/A ff . The graviton mass 
is determined by the ratio k/Mpp Our results are presented in the plane of me, A^ to 
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allow comparisons to be made with any model in this class. We have also defined a specific 


test model (identical to that used in [ 17|), with a value of k/Mpi = 0.01 (at the bottom 
of the range suggested in Q|), which according to |HJ is on the edge of 95% exclusion 
for a first graviton excitation mass of less than 2000 GeV. Thus, we assign a low coupling 
constant to the graviton, and hence obtain a conservative estimate of the production cross 
section. This choice leads naturally to a narrow resonance. This test scenario is used to 
illustrate the potential physics reach in each decay channel. As already emphasised, the 
results derived do not depend on the validity of this particular scenario, but can be applied 
to any model giving rise to narrow well-spaced graviton resonances. For example, our 
results in Section [| show that in a model with k/Mpi = 0.1 (A^ = 10 TeV), the discovery 
limit in the e + e _ channel rises to 3.5 TeV, as a consequence of the increased production 
cross section. Values of k/Mpi > 0.1 are disfavoured on theoretical grounds because the 
bulk curvature becomes too large |p0| . 

An event generator capable of simulating the production and decay of spin-2 resonances 
has been developed. This generator is now part of the standard HERWIG |27], 281 simulation 
package (versions 6.2 and later). The generated events are passed through the ATLAS fast 
simulation (ATLFAST [p9|] ), in order to give a realistic description of detector resolution 
and efficiency. 

In the following sections, the event generator is described (section |2|), followed by 
studies of graviton decays to leptons (section [|), photons and massive vector bosons (section 
|5|), hadronic jets (section |[) , and Higgs bosons (section [/]). Finally, the ability of the LHC to 
determine model parameters, including the length scale of the extra dimension, is discussed. 
Here again we use the simplest RS scenario for illustration. 


2. The event generator 


The implementation of the graviton resonance in the HERWIG event generator has been 
described in [17|. The graviton decays are treated as 2 —> 2 processes, consisting of the 
two hard production subprocesses qq —> G and gg —> G, followed by the graviton decay. 
The relevant matrix elements were computed from the Feynman rules given in (3^, |3l| . 
Interference with Standard Model background processes is neglected. For the range of 
parameter values considered here, the resonance is so narrow that its observed width is 
determined by the detector resolution in all decay channels, and hence interference effects 
cancel in all observable distributions. Note, however, that the neglect of interference is not 
a good approximation for the broad resonances considered in Jlgj ]. 

The production cross-sections in [FTi]] were calculated using the parton distribution 


functions (PDFs) of Owens 132], set 1.1. The present work uses the more recent MRST [33] 


PDFs. This change has no effect on the conclusions of [17]. 

The resulting graviton production cross section at the LHC for the test model is shown 
as a function of the graviton mass in Figure |]. The dashed curve shows the predictions for 


the ‘central gluon’ and ‘high gluon’ leading-order PDFs of [33|; the solid curve is for the 


average of these, which is the default PDF set for HERWIG version 6.3 [34]. This has been 


found to give the best agreement with recent next-to-leading order fits [35|. The dashed 
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curves give an indication of the uncertainty due to reasonable variation of the gluon PDF. 
We see that, even with our conservative choice of the coupling, for a graviton mass of 1.5 
TeV the expected number of produced gravitons is about 5000 for an integrated luminosity 
of 100 fb _1 , falling to about 70 at a mass of 3 TeV. 



Figure 1: Cross section for graviton production at LHC. 

The dot-dashed curve in Figure [T] shows the gluon fusion contribution to the cross 
section for the default PDF set; this is shown as a fraction of the total production cross 
section in Figure Gluon fusion dominates the cross section for graviton masses up to 3.4 
TeV. This has important implications for the angular distribution of the graviton decay 
(see below). 

The branching fractions of the graviton into various decay modes are shown in Figure |3|. 
These predictions are rather model-independent, depending only on the universality of the 
coupling. We see that decays into quark and gluon jets will predominate, due to their 
high multiplicity of colour, spin and flavour states. The Higgs boson fraction depends 
significantly on the assumed Higgs mass when me < 10m#; we have used niH = 115 GeV, 
the default HERWIG value. Out of 5000 produced gravitons with mass 1.5 TeV, we expect 
roughly 3500 jet-jet, 100 e + e _ , 100 100 t + t~, 300 vu, 200 77 , 450 W + W~, 225 

Z°Z° and 15 H°H° decays. 

The angular distributions of the various decay modes in the graviton rest frame are 
summarized in Table |]. Here (3 represents the velocity of the decay products, /3 = 
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Figure 2: Contributions of gluon-gluon and quark-antiquark fusion to graviton production at LHC. 


1 — 4 m 2 /mQ for particles of mass m. In Table [lj, the plot letters refer to Figure |4], 
which shows the distributions in the limit of negligible mass ((3 = 1). Note that the angu¬ 
lar distribution depends strongly on the production mechanism. As we saw above, gluon 
fusion predominates, but the contribution of quark-antiquark fusion has a structure that 
tends to flatten the decay angular distribution. Notice that the angular distributions of 
the massive gauge bosons W and Z are slightly different from those of the massless 7 and 
gluon, even in the limit f3 —> 1, owing to their extra longitudinal polarization state, which 
has the same distribution as the Higgs boson. 


3. Measurements of the graviton couplings 


In 117] it was shown that the graviton resonance can be detected up to a mass of 2080 GeV 
in our test model, using the process pp —► G —> e + e~. The limits are model independent 
as long as the graviton couplings are universal and give rise to narrow resonances, with 
widths less than the experimental resolution. The angular distribution of the lepton pair 
can be used to determine the spin of the intermediate state. In our test model, the angular 
distribution favours a spin-2 hypothesis over a spin-1 hypothesis at 90% confidence for 
graviton masses up to 1720 GeV. In this work, we consider the full range of parameter 
space. In some cases, the search reach can be much higher than in the test model. 
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Figure 3: Graviton branching fractions. 
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Table 1: Angular distributions in graviton production and decay. 9* is the polar angle of the 
outgoing fermion in the graviton rest frame. The letters in the “plot” column refer to the curves in 
Figure [|. 


In the following, it is assumed that the graviton will be detected in the e + e - channel, 
with a significance of greater than 5a. This study then assumes that the graviton mass is 
known, which allows signals with significances as low as 3er to be used in the determination 
of the couplings. 
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Figure 4: Angular distributions in graviton production and decay (/? = 1). 

4. Graviton decays to leptons 

4.1 G -> e + e~ 

This channel offers the best chance of discovery of a graviton resonance at the LHC, 
by virtue of the relatively small background from Drell-Yan processes, and the excellent 
mass resolution provided by the ATLAS electromagnetic calorimeter. More details of the 
proposed discovery search can be found in |l7j]. 

Measurements of the graviton coupling were studied by simulating signals at graviton 
masses between 0.5 and 4.0 TeV. These signals were superimposed on the expected back¬ 
ground from Drell-Yan. Electrons were selected with pt > 5 GeV inside the acceptance 
of the ATLAS tracking detector (|r/| < 2.5), using the standard electron reconstruction 
algorithm of ATLFAST, which accounts for the effect of nearby particles on the calorimeter 
signature. The pair with the highest pt were used to construct the graviton mass. The 
mass distribution of the electron pair is well fitted by a Gaussian signal on a background 
of the form arrive , where m ee is the mass of the electron pair, and a and f3 are free param¬ 
eters. The acceptance of the detector varies from 91 to 76% across the mass range, with 
an estimated systematic error of < 1%, and a negligible statistical error. The efficiency 
for detecting an isolated electron is taken as 90%. The systematic error on this value will 
depend on the details of the detector and reconstruction code used, and is hence beyond 
the scope of this study. However, we note that, for the very high energy electrons involved 
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in these decays, it will not be necessary to make tight cuts on the electron tracks, and so a 
high efficiency with a small error should be obtained. The fit for a graviton mass of 1500 
GeV is shown in Figure |j. 



Figure 5: The number of events per 4 GeV mass bin from a graviton resonance, with tog = 1.5 
TeV (signal), superimposed on the expected Standard Model background (SM), for 100 fb _1 of 
integrated luminosity. The fit to the data is shown by the dotted curve. 

In order to estimate the precision which could be expected for the measurement of 
the production cross-section times branching ratio a.B , a procedure for subtracting the 
background under the peak is required. To avoid assumptions about the background shape, 
we use a simple background subtraction procedure. The background estimate Ng st is 
obtained by counting the number of events in two bins of width w/2 on either side of the 
signal. This procedure will work well since the mass window used to select the signal, 
of width w, is narrow. It contains N events, made up of Ns signal and Nb background 
events. The error on the signal estimate N$ st is then given by 

A N! st = yjAN 2 + ANg st 2 . (4.1) 

The fractional error on a.B is then simply equal to A Ng st /Ng st 

In a real experiment, it would be possible to use a more sophisticated procedure to 
obtain , by fitting over a larger range of electron pair mass. This means that our 
estimates of the experimental reach are conservative. But since, in the interesting regions, 
the background levels are small, the effect of background subtraction on the final error is 
also small. 

































A systematic error arises from the errors on the acceptance, electron identification 
efficiency and luminosity. The method for luminosity measurement in ATLAS is not yet 
decided, but an error of 5-10% is within reach of conventional methods [36]. Methods to 
improve this to 1-2% are under consideration. Together with the as yet unknown electron 
reconstruction efficiency (see above) this means that the systematic error on cross-section 
measurements is very uncertain at present, and hence we plot our results using statistical 
errors only. The effect of the systematic error on the extraction of model parameters is 
discussed in section ||. 

The above procedure was used to determine how well the ATLAS detector could mea¬ 
sure the graviton coupling. The HERWIG/ATLFAST simulation was run at each graviton 
mass, and for the standard model background, for 100 fb~ 4 of integrated luminosity, cor¬ 
responding to one year of running of the LHC at its nominal luminosity of 10 34 cm _2 s _1 . 
The estimated statistical errors on a.B are plotted on Figure ||. Also shown are the re¬ 
gion excluded by Tevatron data [37, 38] and lines of constant k/Mpi. In the test model, 
with k/Mpi = 0.01, a 10% measurement of the production rate is possible for graviton 
resonance masses as high as 1400 GeV. In models with k/Mpi = 0.1 (A*- = 10 TeV), 
a 20% measurement of the coupling is possible for graviton masses as high as 3.5 TeV, 
indicating the ultimate search reach. For completeness, we have continued the contours 
into the region of very small k/Mpp This region is excluded in the RS—I model [^] , but 
may be relevant to other models with low mass graviton resonances. In such a case, values 
of A^ as high as 100 TeV are in principle accessible. For high values of k/Mpi > 0.05, 
horizontally-striped on the plot, the graviton resonance width becomes larger than the 
experimental mass resolution, making a measurement of the width possible. It is unlikely 
that such a measurement would be possible in any of the other channels considered here, 
since their mass resolutions are far inferior. Once the width becomes very large, it would 
be necessary to increase the size of the mass window in order to maintain a high efficiency 
for the signal. Interference with the Drell-Yan background would also need to be taken 
into account. However, the production cross-section is proportional to ( k/Mpi ) 2 , and so 
the detection of the resonance is trivial in this region. 


4.2 G -*• p + p- 

The analysis of this channel is very similar to that for the electron case. Muons were selected 
with pt > 5 GeV using the standard muon reconstruction algorithm of ATLFAST. The pair 
with the highest pp were used to construct the graviton mass. The discovery potential in 
the p + channel is not as great as in the e + e _ case, because the momentum resolution 
of the magnetic spectrometer decreases at high muon momentum. This is reflected in a 
much poorer mass resolution for high muon pair masses. In addition, the reconstruction 
efficiency is poorer than for electrons. Nonetheless, precisions of a.B < 10% are achievable 
in the test model case, for graviton resonance masses up to 1250 GeV. Figure [7] shows the 
results, which can provide a valuable check of lepton universality in the graviton couplings. 


4.3 G t + t 

The t + t~ decay mode would be extremely hard to observe on the large standard model 
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-^-= 0.01 



Figure 6: Contours showing the statistical precision expected for a measurement of a.B in the 
decay mode G —> e + e _ for 100 fb^ 1 of integrated luminosity (Solid lines with markers). Also shown 
are lines of equal k/Mpi (dashed). The test model has k/Adpi = 0.01. The (red) blocked region 
near the origin is excluded by Tevatron data js?], |3^|. In the (green) horizontally-striped region 
above k/Mpi ss 0.05 the resonance width is larger than the experimental resolution. The (blue) 
blocked region with k/Mpi >0.1 is disfavoured theoretically. 


background from QCD jets. The missing energy from the r decays would spoil the mass 
resolution for the graviton signal, further degrading the significance of any peak. This 
signal is therefore not considered further. 

5. Decays to vector bosons 

The detection of decays to vector bosons would be very important in establishing the 
nature of a graviton resonance, since the universal coupling would be very different from 
that expected for other exotic objects, such as a Z'. In addition the angular distribution 
of the decay products is a characteristic signature of the resonance spin (see Table |I] and 
Figure ||). 

5.1 G —> 77 

The method used to study this channel is identical to that used for the electron case. The 
minimum photon pp was set to 1 GeV. A photon detection efficiency of 90% was applied 
after the effect of the standard photon selection cuts in ATLFAST. 
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Figure 7: Contours showing the statistical precision expected for a measurement of a.B in the 
decay mode G —> /x + for 100 fb -1 of integrated luminosity, as for Figure 0. 


The mass resolution for photon pairs is excellent, being very close to that for electrons. 
However the background is much less well understood. The HERWIG calculation of the 
cross-section only includes the Born term and production from gluon-gluon interactions 
via a quark box diagram. It is known that these diagrams alone predict a cross-section 
for photon pair production at the Tevatron which is a factor of ~ 5 too small [39]. This 
large discrepancy means that one cannot rely on existing Monte Carlo simulations to pro¬ 
duce a reliable background estimate, and in particular that the angular distribution of the 


background cannot be trusted. We note that these effects were not considered in [40]. 

As discussed in Section 2 , graviton production at accessible masses is dominated by 
gluon-gluon fusion, and therefore the photon angular distribution is strongly forward- 
backward peaked, as shown by curve c in Figure [|. An important background, not presently 
simulated by HERWIG (or any other current event generator), is bremsstrahlung from ini¬ 
tial state partons, which is also strongly forward-backward peaked. For these reasons, we 
do not attempt an analysis of how well the resonance spin could be determined at the LHC. 
However, it should be noted that the background level and angular distribution can easily 
be measured in the experiment by using data away from the resonance itself. The photon 
channel will then be very significant in establishing the nature of the resonance. 

For this analysis, we use a photon pair background 5 times as large as predicted by 
HERWIG in order to estimate the precision which could be reached on the coupling in this 


- 11 













channel. The results are shown in Figure 




^-= 0.02 

M P | 

^= 0.03 
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-^-=0.05 


Figure 8: Contours showing the statistical precision expected for a measurement of cr.B in the 
decay mode G —> 77 for 100 fb _1 of integrated luminosity, as for Figure ||. 

5.2 G -»• WW -> Ivjj 

Graviton decays to a pair of W bosons are best detected in the mode where one W decays 
leptonically to the electron or muon final states, and one hadronically. This mode has a 
reasonable branching fraction (29% of W pairs), and the reconstruction of the events can 
be performed by assuming that the missing energy is due only to the neutrino from the 
leptonic decay and using the W mass as a constraint. 

This channel has considerable difficulties compared to the leptonic and 77 channels 
due to the large background from tt and W +2 jets. Consequently this channel would not 
be a discovery channel, but would be useful in confirming the universality of the graviton’s 
coupling. Scale factors for the background were obtained by comparing the HERWIG cross- 
section with the NLL cross-section for tt [[II]] and by using the prescription in ^2| for W +2 
jets, giving scale factors of 1.8 and 1.7 respectively. Background from WW production was 
also included for completeness but was negligible compared to the other backgrounds. 

Standard ATLFAST settings for high luminosity were used except for the jet recon¬ 
struction algorithm, for which the Mulguisin algorithm |43|] with a minimum distance of 
A R = 0.2 was used. This algorithm was chosen because it was found to give the highest 
signal reconstruction efficiency. 
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Jet reconstruction is problematic because the jets come from a W which is highly 
boosted. This naturally leads to two jets very close in rj — (f>: for a 3 TeV graviton, both 
jets often fall within the same calorimeter cell as defined by ATLFAST. 

The neutrino 4-momentum is reconstructed from p and fixing the mass of the Iv 
system to be the mass of a W. This gives a quadratic equation which is solved and the 
average of the two solutions taken. The second W is reconstructed from the highest pp 
jet and a jet which gives 65 < Mjj < 85 GeV. If more than one possibility is found, the 
combination that gives the highest pp{W) is chosen. 

The pt cuts are set fairly low to allow the whole range of graviton masses above 500 
GeV to be treated in a single analysis. Some improvement could be made by tuning these 
cuts for a smaller mass region, using the mass of the resonance as measured in the e + e - 
channel as a guide. 

To reduce the it background the following cuts were imposed: a top reconstruction 
veto which attempts to reconstruct a top mass from each of the IFs and an acceptably 
close jet (within A R < 1), and a cut on the number of central (|? 7 | < 2) jets. 

The signal reconstruction efficiency is 22% at 1.5 TeV, dropping to 6% at 3 TeV. 
The mass resolution is 6%. The background is not well described by any simple form 
over the entire range, but it is very smooth. Therefore it is expected that background 
subtraction will be successfully achieved by fitting a function in the sidebands. The error 
on the background estimate was found by fitting an exponential with a window of m ± 2cr 
in a fit region m ± 6<r. Figure [9| shows the background subtraction procedure, in the case 
k/M Pl = 0.05. 


The error on the signal estimate, Ng st is given by equation (4.1), taking the error on 
to be the statistical error on the number of events in the fit region. The estimated 
statistical errors for a.B were then calculated in the same way as the other channels and 


N% st 


are shown in Figure 10 


In the test model with k/Mpi = 0.01, a measurement of a.B is not possible. However, 
for higher values of k/Mpi, measurement is possible with statistical accuracies in the range 
5-30%. For high masses, the overlap of jets means that the efficiency of the selection cuts 
becomes very low. A different analysis would be required, perhaps based on event shapes 
in order to extend the reach further. At low masses, the effects of the pp cuts on the 
determination of the fit parameters can be seen. 


5.3 G —> ZZ —> lljj 

This channel is analyzed in much the same way as the WW channel. The principal differ¬ 
ence is that the dominant background is Z+2 jets - there is no tt equivalent. Consequently, 
similar reconstruction and cuts were employed with the omission of the cuts intended to 
reduce the tt background. 

The same background fitting and subtraction procedure as in the WW case is followed 
leading to a similar error determination, the results of which are shown in Figure Hi- 

Signal reconstruction efficiency is similar to the WW case and mass resolution is better, 
~ 3%. The smaller branching ratio for this decay is offset by the lower background and the 
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Figure 9: A fit to the background in the sidebands around the G —> WW signal (solid line) and 
its extrapolation under the signal peak (dashed line). 

final reach is very similar to the WW case. As before, the effect of decreasing efficiency at 
high mass is evident as is the effect of the pt cuts at lower masses. 

6. Decays to hadronic final states 
6.1 Inclusive decays to 2 jets 

The signature for the decay mode G —> jj comprises two energetic jets in the detector, 
producing a large transverse hadronic energy. The dominant background process to this 
topology is QCD multi-jet production, which forms a continuum irreducible background. 
This large background would make it impossible to find a statistically significant signal 
peak without the knowledge of the peak position provided by other channels such as the 
e + e” final state, except in cases with very large graviton couplings. However, as stated 
above, we assume that any discovery will be made in other channels, and investigate the 
potential to measure the coupling strength alone in hadronic channels. 

G —> jj candidates are selected by requiring at least two jets with minimum transverse 
energies (Ep) of a quarter of the graviton mass. The continuum background is shown in 
Figure |l2|a, with a resonance at 1500 GeV superimposed, using k/Mpi = 0.08. The signal is 
not visible to the eye. The overall acceptance for the signal selection cuts ranges from 40% 
to 60% depending on the model parameters me and A^. Although the signal observability 
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Figure 10: Contours of statistical precision for a.B in the G 


WW channel, as for Figure 


is not sufficient for a discovery in this decay channel, it is adequate for measurement of 
the graviton coupling over much of parameter space; Ng/y/Wp varies from 3.8 (for me = 
500 GeV, k/Mpi = 0.02) to 0.1 (for mg = 2200 GeV, k/Mpi = 0.1). Ng and Np are the 
number of G —► jj signal and background events after the selection cuts. 


The mass of the graviton is determined from the invariant mass of the two highest Ep 
jets; if a third jet is in close proximity (Ar/, A <f> < 1) to one of these two high Ep jets, the 
mass of the graviton is calculated from the invariant mass of the three highest Ep jets. The 
mass resolution for this decay mode degrades dramatically with a long tail toward lower 
masses ( e.g ., the mass resolution is about 160 GeV for a 2000 GeV graviton). 


To extract the signal, a procedure independent of the signal shape has been adopted. 
The mass distribution of the jets is fitted by a Gaussian signal (whose peak is fixed to 
be at the previously determined peak of the signal) on a background of exponential form. 
The procedure for subtracting the background under the peak is the same as that for 
the G —> WW channel. Figure 12b shows the signal after background subtraction. The 
very high statistics in the background sidebands allow the background under the signal 
to be estimated with very high precision, revealing the signal peak, 
fractional error on a.B versus A„- and graviton mass in the G 
of integrated luminosity. 


Figure 13 shows the 
jj decay mode for 100 fb^ 1 
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Figure 11: Contours of statistical precision for a.B in the G — * ZZ channel, as for Figure |j. 

6.2 Decays to heavy quarks 

We have considered whether tagging heavy quarks (t or 6) might be useful to improve the 
signal to background ratio and allow a better measurement of the graviton coupling to 
quarks. However, since both the graviton and gluon couplings to quarks are flavour inde¬ 
pendent, no extra discrimination is obtained. The best experimental strategy is therefore 
to exploit the higher statistical power of the inclusive dijet channel. 

7. Decays to Higgs pairs 

In the case that me > 2m h then the channel G —> HH is open. The angular distribution 
is given in Table [I]. The branching ratio is 1/12 of that into di-photons, if me 3> 2 mu- 
The dominant final state will be 4 6-quarks (for a light Higgs boson) or 4 IT’s (for a heavy 
Higgs boson). Both of these final states are difficult to reconstruct, and will have poor 
statistics. The Higgs coupling is therefore unlikely to be measurable at the LHC. 

8. Summary of the reach for coupling measurements 

The results of the simulations of each graviton decay channel are summarised in Table 
|2[ The precision which can be reached on a.B in each of the channels investigated, is 
presented for a range of points in the me, A*- plane. 
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Figure 12: Distribution of invariant mass a) for all two jet events selected, and b) after background 
subtraction. 


Channel 

Point mo , A n (TeV) 

1,10 

1,20 

1,30 

2,10 

2,20 

2,30 

3,10 

3,20 

e + e _ 

1.6 

3.3 

5.3 

5.4 

11.0 

17.1 

15.1 

30.7 


1.9 

4.5 

8.2 

6.2 

15.2 

28.2 

15.1 

32.7 

77 

1.2 

2.9 

5.2 

3.9 

8.8 

15.2 

10.5 

23.0 

WW 

11.6 

44.9 

- 

38.2 

- 

- 

- 

- 

zz 

13.7 

50.1 

- 

52.7 

- 

- 

- 

- 

33 

19.0 

77.0 

- 

31.0 

- 

- 

59.0 

- 


Table 2: The relative precision achievable (in %) for measurements of a.B in each of the channels 
considered, for fixed points in the tog, plane. Points with errors above 100% are not shown. 


9. Determination of the model parameters 


Several models have been built || 10] in which our analysis applies. They are based on the 


original RS model but with additional branes. Supersymmetric versions [14, 15] have also 
been constructed, in which the graviton resonances are identical to the ones studied here. 
These models all have a narrow relatively strongly coupled first massive graviton mode, 
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Figure 13: The fractional error on a.B versus and graviton mass in G —> jj decay mode for 
100 fb -1 of integrated luminosity, as for Figure |[ 

and our analysis should apply to them (as long as the parameters are such that the graviton 
is narrower than the experimental resolution). This is because graviton modes couple to 
matter in proportionality to the energy-momentum tensor as a model-independent feature, 
guaranteeing the universality of the coupling. The overall coupling strength is a model 
dependent parameter, as of course is the connection between the model parameters and 
the mass or coupling of the first graviton mode. We note that even in the factorisable extra 
dimension case [fj 2 [ ^3| , the resonances would be well separated if the extra dimension were 
small enough. The splitting between the resonances in the factorisable case is constant, 
being 2n/R where R measures the size of the extra dimension(s) [|30|. If 1 /R were in 
the range Mz — 1 TeV, resonance graviton production would still not be possible because 
each state couples with negligible strength, suppressed by the Planck mass. But if models 
were constructed which increase this coupling in the factorisable dimension scenario, our 
analysis would be fully applicable to this case also. 

As a specific example, we have considered the precision which could be obtained on the 
parameters of the RS model, in the case me = 1500 GeV. Since the model only contains 
two parameters (mg and = 39 TeV), only two measurements are required to fully 
constrain the model, me can be measured directly in the e + e~ channel, with a statistical 
precision of better than 1 GeV (our fit gives 0.7 GeV for me = 1500 GeV). The energy 
scale error is given in ]36|] as < 0.7% in this energy range, giving a resolution on me of 
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10.5 GeV. The statistical error on the coupling depends on A^, ranging from 1% at low 
A,r to 15% for A ?r =39 TeV as in the test model. The dominant systematic error on the 
coupling is due to the luminosity measurement, since the systematic errors on efficiencies 
and acceptances will be below the 1% level. We assume, conservatively, that the luminosity 
can be obtained to 10%, giving an overall error on a.B of 18%. We can then infer a value 
for the compactification radius of the extra dimension, r c , and its error, using equations 
Q) and ([□]) . The precision on the coupling measurement is then directly reflected in the 
error in r c giving r c = (82 ± 7) x 10~ 33 m. This reach to extremely small distance scales 
is a consequence of the warp factor in the model, working on the TeV-scale measurements 
of physical observables. 

10. Conclusions 

The LHC detectors will be capable of discovering narrow graviton resonances predicted in a 
range of models with extra space dimensions. Such resonances will most easily be detected 
in the di-electron and di-photon final states. The coupling strength of the resonance to 
p + , W + W~, Z°Z° and jet-jet final states (but not t + t~ or H°H°) can also be measured 
over a wide range of parameter space. The resonance spin can also be measured over a 
more limited mass range. Taken together, these measurements would provide compelling 
evidence for the existence of a massive graviton resonance coupling to the SM fields with 
a universal coupling strength. 

Since models with a small number of Planck-scale extra dimensions are highly con¬ 
strained, the model parameters can be extracted with good precision, if a particular sce¬ 
nario is assumed. In the RS test model, the size of the extra dimension can be inferred to 
better than 10%, corresponding to a precision in length of 7 x 10” 33 m, using measurements 
of a.B and the graviton mass. 
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